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Abstract 

In quantum coding theory, stabilizer codes are probably the most important class 
of quantum codes. They are regarded as the quantum analogue of the classical linear 
codes and the properties of stabilizer codes have been carefully studied in the literature. 
In this paper, a new but simple construction of stabilizer codes is proposed based on 
syndrome assignment by classical parity-check matrices. This method reduces the 
construction of quantum stabilizer codes to the construction of classical parity-check 
matrices that satisfy a specific commutative condition. The quantum stabilizer codes 
from this construction have a larger set of correctable error operators than expected. 
Its (asymptotic) coding efficiency is comparable to that of CSS codes. A class of 
quantum Reed-Muller codes is constructed, which have a larger set of correctable error 
operators than that of the quantum Reed-Muller codes developed previously in the 
literature. Quantum stabilizer codes inspired by classical quadratic residue codes are 
also constructed and some of which are optimal in terms of their coding parameters. 



Index terms: Quantum error- correcting codes, quantum stabilizer codes, quantum 
information theory, Reed-Muller codes, quadratic-residue codes. 



1 Introduction 



The theory of quantum error correction has been profoundly developed in the last decade 
since the first quantum error-correcting code proposed by Shor ^J. In [2], using a different 
approach from that of Shor, Steane gave a new quantum error-correcting code and studied 
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basic theory of quantum error correction. Later, Steane gave more new quantum codes and 
discussed constructions of quantum error-correcting codes in [3j. The complete quantum 
error correction condition was given in [UEIE]- The optimal five qubit code was discovered 
in [5l [7]. After CSS code construction [H |9], the study of quantum error- correcting codes 
then turned to the study of classical self-orthogonal codes. 

The idea of stabilizer codes was proposed in [IHIIII] and the properties of stabilizer codes 
were extensively addressed in CSS codes can then be viewed as one prominent class 

of stabilizer codes. In p^, Steane gave a further improvement, called an enlargement of CSS 
codes, which produces several families of quantum codes with greater minimum distance. In 
this paper we will propose a new but simple construction of quantum stabilizer codes based 
on syndrome assignment by classical parity-check matrices and develop several classes of 
quantum codes from this construction. 

This paper is organized as follows. In Section [21 we begin in describing the basic proper- 
ties of stabilizer codes and end with a new formulation of CSS codes and their enlargement 
[El El [13]. The method of construction of stabilizer codes based on syndrome assignment 
by classical parity-check matrices is proposed in Section [3l where we discuss the asymptotic 
coding efficiency of this method. As an illustration, we develop a family of quantum sta- 
bilizer codes from classical Reed-MuUer codes in Section [H Several other quantum codes 
inspired by classical quadratic residue codes are investigated in Section [5l where several 
optimal quantum codes are constructed. A conclusion is discussed in the last section. 



2 Stabilizer Codes 

2.1 Stabilizer Groups and Stabilizer Codes 

Let Ti be the state space of a qubit. The Pauli group acting on the state space 7i®" of n 
qubits, plays an important role in the construction of n-qubit stabilizer codes. An element 
in Qn is expressed as i'^Mi ® M2 ® . . . M„, where each Mj is one of the Pauli operators /, 
X, F, or Z on 7^, i = and c G {0, 1, 2, 3}. Let K, = {±7®'', ±i/®"}, which is a normal 

subgroup of Qn and will be used in a later discussion. 

For a. g E Qn, the fixed subspace V{g) of (7 is a subspace of 7i®" such that \ip) G V{g) 
if and only if 5']^') = A stabilizer group S that fixes a non-trivial subspace T of Ti®^ is 
the set 

S = {gegn\gm = \i^). V|^)Gr}. 
A necessary condition is that — / ^ S. Since any g,h & Qn have either gh = hg or gh = —hg, 
S must be an abelian subgroup of Qn- And since any g E Gn has either g"^ = I 01 g"^ = —J, 
we have g'^ = I Wg & S. Therefore, S = {Jj-if for some r, i.e., S =< gi,g2, ■ ■ ■ ,gr > with r 
commutative independent generators. 

An [[n, k, d]] quantum stabilizer code C{S) is a 2'^-dimensional subspace of 7^®" fixed 
by a stabilizer group S with a set of r = n — independent generators. The d means the 
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minimum distance of the quantum code C{S) and will be defined later. The error- correct ion 
condition for a stabilizer code C{S) [HI [T^l [H] says that {Ei} is a collection of correctable 
error operators in Qn for C{S) if and only if 

E]Ek^N{S)\SWj,k, (1) 

where 

N{S) = {ge Gn \9hg^ eSyheS} 
is the normalizer group of S in Qn, which in fact is the centralizer group of S in and 

S = S}C = {gh \ g eS, h e IC} . 

The weight of an element i'^Mi®M2 ® . . . ® M„ in is defined to be the number of Mj's not 
equal to /. Then the minimum distance d of C(5), motivated by the above error correction 
condition, is defined to be the minimum weight of an element in N{S) \ S. 

2.2 Binary Codes Corresponding to Stabilizer Groups 

The Pauli group Qn is closely related to the 2n-dimensional binary vector space Z^". If u,v 
are two binary n-tuples, (m, v) is meant to be a binary 2n-tuples and any element x G Z2" 
can be written in the form (m, v) with m, v n-tuples. We use uv to denote the n-tuple of the 
bitwise AND of u and v. That is, {uv)i = Ui ■ Vi where the subscript i means the i-th bit of 
the binary ra-tuple. Then we define the generalized weight of an 2?7,-tuple x = {u,v) in Z^", 
denoted by gw{x), as the Hamming weight of the bitwise OR of u and v. Thus 

gw{x) = w{u) + w{y) — w{uv), 

where w{u) means the Hamming weight of u, the number of nonzero components of u. 

A g = i'^Mi ® M2 ® . . . ® M„ in Qn, can be expressed as g = i^' X^Zp, where a = 
(ai, 02, ... , a„) and /3 = (61, 62, • • • , are two binary n-tuples and c, c' G {0, 1, 2, 3}. In this 
expression, if Mj = I,X,Z,Y, then {aj,bj) = (0, 0), (1, 0), (0, 1), (1, 1), respectively. And 
we have c' = c + I (mod 4) where / is the number of Mj's which are equal to Y (note that 
Y = iXZ). We define a group homomorphism r : Gn/^ by 

^{9) = 9^- 

li g = i'^XaZ^, T{g) = X^Zp K. Note that X^ZplC = Xa'Z^ilC if and only if a = a' and 
(3 = (3' . Also r is an epimorphism. Next we define a group isomorphism : Qnl^ ^ ^2" tiy 

/i(X„Z^/C) ^ («,/3). 

Then we can define a homomorphism if : Z2" by = /i o r, i.e., 

i^ii'X^Zp) ^ fi{T{t'X^Zf,)) = /x(X,Z^/C) = {a, (3). 

It is clear that ip is an epimorphism with kernel K. and {5^1, (72, • • • , S'r} is a set of independent 
generators if and only if ip{gi) , {p{g2) , . . . , (p{gr) are linearly independent 2n-tuples in Z^". 
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Suppose {gi,...,gr} is a set of independent generators of a stabilizer group S. We 
define a check matrix if of iS by making (p{gi) as its i-th row vector. Then if is an r x 2n 
binary matrix. For convenience, we may denote H by 



H=[Hx\Hz], 

where Hx, Hz are two r x n binary matrices. Note that each of the 2'' sets {±gi, . . . , ^gr} 
can be used as a set of independent generators of a distinct stabihzer group. But all of them 
have the same check matrix H. 

One example of stabilizer codes is a five-qubit code. A set of generators of a [[5, 1,3]] 
code can be 

^1 = XZZXI, 
g2 = IXZZX, 
gs = XIX ZZ, 
^4 = ZXIXZ. 

The corresponding check matrix is 

"10010 01100" 

01001 00110 , , 

10100 00011 ■ ^ ^ 

01010 10001 



Since a stabilizer group S is an abelian group, we have gh = hg \/g, h E S, which has a 
corresponding property in ip{S), induced by the homomorphism if, that 



where A 



2n 



Or 



Or, 



(p{g)K2nV{h) = 0, V^f, he S, 



. Thus a check matrix H of a stabilizer group S has to satisfy 



the following commutative condition, 



HxHz + HzHx 



Or 



(3) 



where Oixj is the i x j zero matrix. We will omit the subscripts of A and O in the following 
discussion. We define that an r x 2n binary matrix H = [Hx \Hz] is commutative if it 
satisfies the commutative condition. From ([3]), an r x 2n binary matrix H = [Hx \Hz] is 
commutative if and only if HxH^ is a symmetric r x r matrix. 

The check matrix is a convenient tool for the encoding and decoding of stabilizer codes 
jl2] . In addition, the check matrix is able to facilitate the construction of stabilizer groups 
from known classical binary codes, as will be demonstrated in the next section. Before that, 
we will illustrate an application of the check matrix of a stabilizer group to determine the 
minimum distance of the corresponding stabilizer code as follows. 
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Let S = (p{S), which is a subspace of • Then S can be viewed as a classical binary 



2n 



code with generator matrix H. As in |1H IT5]. we define a symplectic inner product * on Z; 

by 

f i) * (U2, V2) = Ui ■ V2 + Vi ■ U2. 

Thus, two elements (7, /i in ^„ is commutative if and only if the symplectic inner product 
ip{g) * (p{h) of ip{g) and (p{h) is zero. Let S-^* be the dual code of S with respect to the 
symplectic inner product, 



{iu,v) G Zf Uu,v) * {a, (3) = 0, V(a,/3) eS}. 



It is clear that 5 is a self-orthogonal code with respect to this symplectic inner product, i.e., 
S C 5"*"*. An (2n — r) x 2n binary matrix 

G = [ Gx I Gz ] 

of rank 2n — r, where and Gz are (2n — r) x n binary matrices, is a generator matrix of 
5*-^* if and only if 



HAG^ 



HxG^z 



HzGx 



O. 



(4) 



It can be seen that S-^* = (p{N{S)) and the minimum distance of C{S) is just the minimum 
generalized weight of a nonzero codeword in S-^* \S. This helps decide the minimum distance 
of a stabilizer code, as illustrated in the construction of CSS codes [8l[9] and their enlargement 
as follows. 



To construct an [[n, CSS code C{S), we choose a classical [n,ki] binary code Ci 
and an [n, ^2] subcode C2 of Ci such that both the code Ci and the classical dual code of C2 
have minimum distance ^ d. Then a check matrix of a stabilizer group S is established as 



' G2 


' 





Hi . 



H 



of rank n — ki + k2, where G2 is a generator matrix of C2 and Hi is a parity-check matrix of 



Ci. Let Gi 



G2 
G3 



a generator matrix of Ci, and H2 



€2- Then a generator matrix of the symplectic dual code S- 



Hi 
* is 



, a parity-check matrix of 



G2 


' 















' Gi 


' 





Hi 







H2 





H3 









G 



which is of rank n + ki — k2. It can be verified that both ([3]) and (jll) hold. The minimum 
distance of the quantum code is no less than the minimum generalized weight of the sym- 
plectic dual code S-^* and is clearly ^ d from the structure of its generator matrix G. The 
dimension of the quantum code is k = n — {n — ki + k2) = ki — k2. 
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The enlargement of CSS codes in [13j is based on the CSS construction and exchanges 
code dimension for error-correcting capabihty. With additional stabilizer generators, this en- 
largement increases the minimum distance of the code by half. We give an explicit description 
of the enlargement in terms of nonsingular matrices and a proof that this enlargement is able 
to increase the minimum distance of the code by half in the following theorem. 

Theorem 1. Let Ci be a classical [n, ki, di] binary code which contains its classical dual C^. 
Furthermore, let Ci be able to be enlarged to C2 = [n,k2,d2], where k2 > ki. Suppose that 

Gil G2 = ^ generator matrices of Ci, C2, respectively, and H2, Hi = ^ are 

parity-check matrices of C2, Ci, respectively. If there exists a (^2 — ki) x {k2 — ki) nonsingular 
binary matrix P such that / + P is also nonsingular, by taking 



H 



H2 








H2 




H3 



where Q = [H^GJ) {P^) ^ [H3GI) ^ is a {k2 — ki) x (/c2 — ki) nonsingular matrix, as a 
check matrix of a stabilizer group S, an [[n, k2 + ki — n, d ^ minjcii, f^]}]] quantum code 
C{S) can be constructed. A generator matrix of 5^"*"* can be 



G 



1 Gi 










PGs J 



Proof. It is easy to see that the two matrices H and G are of full rank and have rank 
{2n — ki — k2) and rank ki + k2, respectively. Since G Ci, we have HiHj = O and hence 
H2HI = O, H3HI = O, H2HI = O. Thus 



HAH^ ={ O \{0^ Hi Hj)+\ H2 \{ Hi HjQ^ ) 




OHlQ^ 

H2hIq^ 

HsHJQ^ 



and Eq. holds. Since HiGj = O and H2G^ = O, we have H2GI = O, H3GI = O, and 
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H2G^ = 0. Thus 



HAG 




GJ 



GlP^ 



H^G^P^ 



OGlP^ 



H2GJ 
OGJ 

QH.GJ QH.GlP^ 
O 

o 




OGl 





( Gi 




codeword (ci,C2) = (mi ^2 M3) 




Gi 






PGs J 



O O {H,Gj){P^)-\H,Gjy'H,GjP^ + H,Gl 



and Eq. (jlj) holds. The minimum distance d of C{S) is determined as follows. Consider a 



{uiGi + u^G^,U2Gi + u^iPG^) in 



S , with two /ci-tuples ui,U2 and a (A;2 — A;i)-tuple U3, not all zero tuples. It is clear that 
Ci, C2 G C2 and w(ci) ^ ^2, w(c2) ^ (^2- If Ci 7^ C2, then w{ci + C2) ^ (i2 and we have 

fifw(Ci, C2) = w{ci) + ti?(Ci) - w{CiC2) 

= ^ + W{C2) + (w(Ci) + W(C2) - 2w{CiC2))) 

= 2 (^('^1) + ^('^2) + w{ci + C2)) ^ — . 

Now if ci = C2, then uiGi + M3G3 = U2G1 + u^PG^ and then (mi + ^2)^1 + M3G3 = U3PG3. 
Since Gi and G3 have linearly independent rows, we must have Ui + U2 = and M3G3 = 
U3PG3. However, since / + P is nonsingular, U3G3 = U3PG3 if and only if U3 = 0. Since 
ui, U2, U3 are not all zero tuples, we have nonzero ui = U2 and ci = uiGi G Ci, C2 = U2G1 = 
uiGi = ci. Thus gw{ci,C2) = w{ci) ^ di. In conclusion, the minimum distance of the 
quantum code is 

3d2 



d ^ min{(ii. 



}• 



□ 



3 A Simple Construction of Stabilizer Codes 
3.1 Syndrome Assignment and Check Matrices 



Given a stabilizer group S together with a set of independent generators {gi, g2, . . . , gr}, 
the encoding-decoding techniques of the corresponding stabilizer code is well studied in 
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But one may ask the question: how to find a set of independent generators gi, g2, . . . , Qr such 
that the generated stabihzer group will correspond to a good quantum code? 

In [T0|, Gottesman constructed a class of stabilizer codes saturating the quantum Ham- 
ming bound which says that 




These codes encode k = n — j — 2'm.n = 2^ qubits and correct up to t = 1 error. He assigned 
an "error syndrome" to each correctable error operator according to a certain rule so that 
error syndromes form a check matrix of the stabilizer group. This is a brilliant idea of giving 
a set of independent generators gi, g2, . . . , gr- We will basically follow this idea to construct 
[[n, k, d]] quantum codes for the general case of t = [^^J > 1. 



3.1.1 Error Syndromes 



In classical coding theory, correctable error patterns are a collection of error patterns which 
have distinct error syndromes. The same concept holds in non-degenerate quantum codes 
and hence can be applied to construct non-degenerate quantum stabilizer codes. In stabilizer 
codes, the idea of error syndrome comes from the commutativity. The error- correct ion 
condition for stabilizer codes in ([T]) says that the multiplication of any two correctable error 
operators in Qn, each with weight no more than t, is not in N{S) — S. So each correctable 
error operator is assigned with a binary pattern indicating the commutative relation between 
the error operator and each one of the generators of a stabilizer code. For our purposes, we 
will take a stricter error- correct ion condition that multiplication of any two correctable error 
operators, each with weight less than or equal to t, must anti-commute with some element 
in S. 



As in [lOj, for any (? G we define : Qn ^ ^2 by 



0, if [g,h]=0, 

1, if {g,h}=0, 

where [g, h] = gh — hg and {g, h} = gh + hg. It can be verified that fg is a group homomor- 
phism. Then for a given set of independent generators gi, g2, gr of a stabilizer group S, 
we define : Qn ^ (^2)'' by 



It is clear that 



fsih) = {Uh),Uh),...,Uih)f 



fs{h) = HA^ihf, (5) 



where H is the check matrix of S corresponding to the (ordered) generators gi, g2, . . . , g^ 
Note that fs is a group homomorphism. It is also clear that fs{h) = (0, 0, ... , 0)"^ if and only 
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if h commutes with every element in S, i.e., h G N{S). We call fs{E) the error syndrome 
of E for each error operator E in Qn- The error syndrome in this definition is equivalent 
to that defined in [13] with +1,-1 in place of 0, 1, respectively. For any two correctable 
error operators Ei and E2 in Qn, each with weight less than or equal to t, we attempt to 
construct a stabilizer group S such that fs{EiE2) is a nonzero vector. Since fs is a group 
homomorphism, 

fs{E,E2)^0^ fsiE,)^ fsiE^). 
As a result, we need to assign distinct error syndromes to distinct correctable error operators. 

If we can correct error operators Xi, . . . , X„, Yi, . . . , y„, Zi, . . . , Z„ and the multiplica- 
tion of any no more than t of them, we are able to correct all error operators up to weight t. 
We call these 3n error operators to be basic correctable error operators. Moreover, since the 
set {Xi, X2, . . . , Xn, Zi, Z2, . . . , Zn} generates all the error operators in Qn by multiplica- 
tion up to a scalar factor, it is desirable to determine the error syndromes of these 2n basic 
correctable error operators so that the error syndromes of all correctable error operators are 
distinct from each other. For the case of t = 1, it has been done in [10]. For the general case 
of t ^ 2, it becomes challenging to assign error syndromes to the 2n basic correctable error 
operators so that the error syndromes of correctable error operators of weight no more than 
t are all distinct. 

From ([5]), we observe that the first n columns and the last n columns of the check matrix 
of a stabilizer group S corresponding to the (ordered) independent generators gi,g2, ■ ■ ■ , dr 
are /^(^i), fs{Z2), ■ ■ -JsiZn) and fs{Xi), fs{X2), ■ ■ •, fs{Xn), respectively. Thus we can 
establish a check matrix of a target stabilizer group S by assigning 2n error syndromes as its 
columns fs{Zi), fs{Z2), ...,fs{Zn), fsiXi), fs{X2), . . ., fs{Xn) and verifying this matrix 
to be commutative. In this way, the method of syndrome assignment is just to define a legal 
check matrix. 



3.1.2 Syndrome Assignment by a Binary Parity-Check Matrix 

Let E = XaZj3 be an error operator of weight no more than t*, where a = (ai, . . . ,a„), 
P = {bi, . . . ,bn) and gw{a, j3) < t*. Then for a target stabilizer group S, 

n n 

fs{E) = fsiX^Zp) = 5^a,/5(X,) + Y^bJs{Zi), 

i=l i=l 

which is a linear combination of at most 2t* terms. To ensure that the error syndromes 
of two distinct error operators Ei = X^^Zj^-^ and E2 = X^^Zp^ with ai = {ai, . . . ,an), 
Pi = (bi, . . . gw{ai,Pi) < t* and a2 = . . .,Un),P2 = {vi, . . .,Vn), gw{a2,P2) < t* 
are distinct, we must have 

fs{Ei) ^ fs{E2) ^ fsm + ^ 

n n 

^ ^(a, + u,)fs{X,) + Y,ib^ + v,)fs{Z,) ^ 0, (6) 
1=1 1=1 
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which is a hnear combination of at most At* terms. Thus a sufficient condition to guarantee 
the distinction among error syndromes of error operators of weight no more than t* is that 
any At* elements in the set {fsiXi), /^(Xa), . . . , fsiX^), fs{Zi), fs{Z2), fsiZn)} must 
be hnearly independent. Surprisingly, this is just a property of a parity-check matrix of a 
classical linear block code with minimum distance d' ^ At* + 1, where any d' — 1 column 
vectors of the parity-check matrix must be linearly independent. 

From the above discussion, we now know how to do syndrome assignment so that all 
correctable error operators will have distinct error syndromes from each other. We ffist 
choose a classical [2n, n + k, d'] binary linear block code C with d' > At* + 1 and k > 0. Let 
H' be a parity-check matrix of C with dimension {n — k) x2n . Then the 2n column vectors 
of H' will be assigned as /5(Xj)'s and fs{Zi)^s. There are (n — k) independent generators 
of the target stabilizer group S and the corresponding target check matrix H is just the 
permutation of the column vectors of H', i.e., H = H'P for some permutation matrix P. 
Then the commutative condition becomes 

H'PAP^H'^ = 0. (7) 

If G' is a generator matrix of C, then the symplectic dual S^* of S in has a generator 
matrix 

G = G'PA 

since 

HAG^ = H'PAA^P^G'^ = H'G'^ = O. 

If the target check matrix H satisfies ([7]), then H is indeed a check matrix of a stabi- 
lizer group which corresponds to a quantum code that is at least t*-error correcting. The 
dimension of the quantum code is n — {n — k) = k. Thus the choice of k > ensures that 
the corresponding quantum code is of dimension greater than zero. 

We conclude the above discussion in the following theorem. 

Theorem 2. Given an {n — k) x 2n parity-check matrix H' of a binary [2n, n + k, d'] linear 
block code C with minimum distance d' ^ At* + 1, such that ([7j) holds for a certain per- 
mutation P, an [[n, k,d > 2t* + 1]] stabilizer code with a check matrix H = H'P can be 
constructed. 

A t*-error-correcting quantum code of length n has 

1=0 ^ ^ 

error syndromes. One may expect that a quantum code by the above construction can 
correct more than just those error operators of weighty t* when t* is determined from d^ 
in Theorem [21 In fact, any error operator E = XaZp with w{q) + w{l3) ^ 2t* has its own 
unique syndrome and then can be corrected. For example, error operators and with 
w{a) = 2t* are correctable. 
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t + 1 








t* + lt + l 



2t* 



mx,z 



Figure 1: The mx,z — region of additional correctable error operators. The mx,z-axis 
represents the number of components M^-'s of error operator E equal to X or Z and the 
my-axis represents the number of components Mj 's equal to Y . 

Let E = i^Mi (g) Ma ® . . . ® M„ = X^Zp, a = (ai, . . . , a„), /3 = (61, . . . , 6„) e , be a 
correctable error operator with w{a) + ^ 2t*. Let / = gw{a,P), the weight of E, and 
then ^ / ^ 2t*. Let my = w{a(3), the number of M,'s equal to Y, and then mx,z = I — my 
represents the number of the M/s equal to X or Z. Since w{a) + w{[3) = 2my + mx,z = 
my + I, we have my + / ^ 2t* (note that I = mx,z + fny)- There are 



correctable error operators for a certain / and a certain my satisfying ^ my ^ / and 
^ my + 1 ^ 2t*. Summing / from t* + 1 to 2t* and summing my from to 2t* — /, we have 
additional 



correctable error operators of weight > t*, which can be a large amount! Figure [T] illustrates 
the mx,z—^Y region of all additional correctable error operators of weight / = mx,z + my > 
t* as the dashed triangle. 

On the other hand, for a given stabilizer group S with quantum error-correcting capa- 
bility t, the classical minimum distance d! of a check matrix of S can be used to determine 
the existence of additional correctable error operators of weight > t and how many of them 
as stated in the following theorem and illustrated in the dashed triangle in Figure [H 

Theorem 3. Let t be the quantum error correcting capability of a stabilizer code C{S). Let 
t* be an estimate of t by a check matrix of S as in Theorem H If t < 2t* or [^J < 2[^J, 
then we have additional 





(8) 
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correctable error operators of weight > t. 



Finally we give a slight improvement of Theorem [2j Define a matrix Hy = Hx + Hz for 
a check matrix H = [ Hx \ Hz ) of S. Let Cx, Cz, Cy be classical binary linear block codes 
with parity-check matrices Hx, Hz, Hy, respectively. If error syndromes of all correctable 
error operators of weight no more than t* are distinct, it is necessary that Cx,Cz,Cy have 
minimum distance ^ 2t* + 1. 

Corollary 4. If Cy has minimum distance ^ 2t* + 1, the condition in Theorem [2] can be 
reduced to d' > At*. 

Proof. We only need to consider two distinct error operators Ei = X^^Zp^ and E2 = X^^Zp^ 
with «! = (ai,...,a„), j3i = (61,..., gw{ai,pi) = t* and 03 = P2 = 

{vi, . . . , Vn), gw{a2, P2) = t* such that 

n n 

fs{Ei) + fsiE2) = ^(a, + u,)fsiX,) + ^(6, + v,)fsiZ,) 

i=l i=l 

is a linear combination of 4t* columns of H. Then we must have ai = j3\, a2 = P2, 
w{ai + 02) = 2t* and then 

n 

fs{Ei) + fs{E2) = $^(a, + M,)/5(y.), 

i=l 

which is a linear combination of 2t* columns of Hy. Since Cy has minimum distance ^ 2t* + l, 
we have 

fsm + fs{E2) ^ 0. 

Thus ([2]) holds for the extreme case — a linear combination of exact At* terms. □ 



3.2 Constructions of Check Matrices 



Our construction of a check matrix of a stabilizer group needs a binary commutative parity- 
check matrix of even length. We suggest three ways to establish commutative parity-check 
matrices by using classical constructions of new codes from old ones [16] such that the 
minimum distances of the resulted quantum codes can be determined from the corresponding 
classical binary linear block codes. 



Construction I: Let Ci, G2 be generator matrices of an [n, fci, di\ and an [n, /c2, 1^2] binary 
linear block codes Ci,C2, respectively, such that ki + k2 > n. Let Hi,H2 be parity-check 



matrices of Ci,C2, respectively. Then C 





' 





C2 



is a generator matrix of a [2n, ki + 



k2, d' = minjc?!, ^2}] code with a parity-check matrix H' 



' H, 


' 





H2 



. H' is commutative 



if and only if H1H2 = O, i.e., the classical dual code of C2 is a subcode of Ci. In this way. 
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' 




G2 






by choosing H = H' as a. check matrix of a stabihzer group S and G = G'A 

as a generator matrix of the symplectic dual code S^*, C{S) is an [[n, ki + k2 — n, d > d* = 
min{di,d2}]] quantum code, where d* equals to the minimum generalized weight of S-^*, 
which corresponds to a vector of minimum weight in Ci or in C2. Since d* = d', we have 
L^J > 2L^J and there is no additional correctable error operators of weight > 
guaranteed by Theorem [31 Note that the construction of CSS codes is a special case of 
Construction I. 

Construction II: Let Gi,G2 be generator matrices of an [n,k,di\ and an [n, A;, ^2] binary 
linear block codes Ci,C2, respectively. Let Hi, H2 be parity-check matrices of Ci,C2, respec- 
tively. Then G' = [Gi \ G2] is a generator matrix of a [2n, k, d' ^ di + ^2] binary linear block 

'Hi O 



where A, B are two matrices such 
lich is greater than n. Since H' has 



code C . A parity-check matrix of C is H' = O H2 

A B 

that GiA^ + ^2^^ = O and the rank of H' is 2n - k, wl 
too many rows to be a check matrix, we consider the code C'-*", the classical dual code of C, 
instead. If the matrix G' is commutative, we choose H = G' = [Gi \ G2] as a check matrix 
of a stabilizer group S. However, it remains to determine the classical minimum distance of 

Construction III: {\u\u + v\ construction) Let Gi and G2 be generator matrices of an 



[n, ki, di] and an [n, k2, ^2] binary linear block codes Ci, C2, respectively. Then G' 



Gi 
O 



is a generator matrix of a [2n,ki + k2,d' ^ minjc/i, ^2}] code. Let Hi,H2 be parity-check 



Gi 
G2 



matrices of €1,62, respectively. A parity-check matrix is H' 



' H2 


H2 ' 


. Hi 






H' is commuta- 



tive if and only if HiHj = O, i.e., the classical dual code of C2 is a subcode of Ci. In this 
case, we take H' as a check matrix if of a stabilizer group S. The minimum generalized 
weight of a generator matrix G = G'A of S-^* is minjdi, ^2}. Hence the stabilizer code C{S) 
has parameters [[n, ki + k2 — n,d ^ mm{di, ^2}]]- 



When Ci is a subcode of C2, we consider the effect of a permutation matrix P' 



I O 
O P 



on the + construction. Let H" = H'P' 



H2 H2P 
Hi O 



H" is commutative if and 



only if H2PH2 = H2P^H2 and H2PHJ = O. In this case, we take H" as a check matrix 



H of a stabilizer group S. Then G = G'P'A 
Consider a nonzero codeword 



GiP Gi 
G2P O 



is a generator matrix of 5* 



{Cl,C2 



(Mi|m2) 



G2P 

GiP 



G2 
O 



{U1G2P + U2GiP,UiG2] 



in S^*, with a fci-tuple Ui and a fc2-tuple U2, not both zero tuples. It is clear that Ci G C2P, 
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C2 G C2 and w{ci) ^ d2, w{c2) ^ ^2- Then 



gw{ci, C2) = w(ci) + w{ci) - w(ciC2) 

= 7; (^(Cl) + W^(C2) + {w{ci) + W{C2) - 2w(ciC2))) 



>d2 + 



'w{ci) + W{C2) + w{ci + C2)) 
w(ci + C2) 



If 



min -w(ci + C2) > 0, 

ci^O or C2^0 2 



we can obtain a quantum code of a greater minimum distance. 

Theorem 5. Let Ci be a subcode of C2 with (ii > ^2 and let ifi, if2 be parity-check matrices 
of Ci,C2, respectively. Assume that P be a permutation matrix such that 



O 



(9) 



and 



Let Gi and G2 



G3 



CiP = Ci, C2P = C2. (10) 
be generator matrices of Ci and C2, respectively, such that 
uG^P 7^ uG^ for any nonzero (^2 — fci)-tuple u. (11) 



Then there is an [[n, k2+ki—n, d ^ minjc/i, [^] }]] stabilizer code C{S) with a corresponding 
stabihzer group S such that 



G = 

is a generator matrix of S-^*. 
Proof. A check matrix of S is 



G2P 
GiP 



G2 
O 



G3P 


G3 


GiP 




GiP 






Hi HiP 



Ho 



O 



since Eq. (|3]) holds by ([9]). Consider a nonzero codeword 



(C1,C2) = (mi|m2|m3) 



/ G3P 


G3 \ 


GiP 




V GiP 


J 



(M1G3P + (m2 + M3)GiP, U1G3 + M2G'i) 
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in S~^*, with two A;i-tuples ui, U2 and a {k2 — A;i)-tuple M3, not all zero tuples. It is clear that 
ci e C2P = C2, C2 e C2 and w{ci) ^ ^2, ^(02) ^ d2. If Mi 7^ 0, 

Cl + C2 = U1G3P + {U2 + U^)GiP + MiGs + 1x2(^1 

= (M1G3P + M1G3) + {{U2 + M3)GiP + M2G'i). 

Since UiG^iP + MiG'3 G C2 - Ci by ([TOD and ^ and (u2 + u^)GiP + ^2^1 £ Ci by (fTO|). we 
have Cl + C2 7^ in C2 and tu(ci + C2) ^ ^2- Then 

1 3(i 

gw{ci, C2) = - (w(ci) + w{c2) + w(ci + C2)) ^ 

If Ml = 0, we have Ci = {u2 + u^jGiP G Ci and C2 = M2G'i G Ci. Hence 

5fw(ci, C2) = w{ci) + w{c2) - w{ciC2) ^ di. 
Therefore, the minimum distance of the quantum code is ^ mm{di, f^]}. □ 

Further investigation on the classical minimum distance may guarantee additional cor- 
rectable error operators of weight greater than the quantum error- correcting capability by 
Theorem [31 A family of quantum Reed-MuUer codes will be constructed by the |m|m -|- 
construction in Section HI 



3.3 Existence of Commutative Parity- Check Matrices 

There is an important question: for a given r x (2n) parity-check matrix H, where r < n, 
does there exist an effective permutation matrix P such that HP is commutative? 

To answer this question, we run a simulation on a computer as follows. Let H = 
[ /rxr I -B ] , where S is a randomly generated r x (2n — r) matrix. Each element of B is 
1 or with probability pi and po = ^ ~ Pi, respectively. By exhaustive search with (2n)! 
permutations for the case n = 5, unfortunately, we found that there exists an H which has 
no effective permutation matrix P such that HP is commutative. However, there is a high 
probability that for a randomly generated matrix H = [ J^xr | ] > there is an effective 
permutation matrix P such that HP is commutative. Moreover, if pi < pq, the probability 
becomes higher. This simulation suggests that parity-check matrices of classical LDPC codes 
may be transformed into legal check matrices by Theorem [2l However, it becomes extremely 
harder to verify this suggestion for n > 8 due to prohibitive computing complexity. The 
question of determining an effective permutation matrix for a parity-check matrix remains 
open. The construction of quantum stabilizer codes can be converted to the construction of 
classical linear codes with commutative parity-check matrices. 

3.4 Asymptotic Coding Efficiency 

In this subsection, we will investigate the asymptotic coding efficiency of the construction of 
stabilizer codes as stated in Theorem[2]by assuming that among all [n, k, d] binary linear block 
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codes, there is at least one code with a parity-check matrix H and an effective permutation 
matrix P such that HP is commutative. 

Suppose an [[n, k,d > d* = 2t* + 1]] stabihzer code with a check matrix is constructed 
by Theorem [2] from a certain [n' = 2n, k' = k + n,d' > At* + 1] classical linear block code, 
where t* = L^^J- Let a' = limsup^/^o^ TP' a = limsup^^^^ K Since 

, , k+n 1 1 

a = limsup = — I — a, 

2n 2 2 

we have 

a = 2a'-l. (12) 
Let 5' = ^,6* = ^ and 5=^. Since I < ^ < I + 1, we have t* <^ <t* + l 

n'' n n L4J — 4 — L4J'' — 4 — 

and then 2d* - 1< d' < 2d* + 3. Thus we have < 6' < and then 6* - ^ < 6' < 

— — 2n — — 2n 2n — — 

(5* + Thus for sufficiently large n, we have 

5* - 5'. (13) 

It is obvious that 

S>5* = 5'. (14) 
From |17|, the classical Hamming bound says that 

a'{6')<l-H,{^-6'), 

where H2{x) = — xlog2(a;) — (1 — a;) log2(l — x). By f|T2|) and f|T3|) . we have a corresponding 
quantum Hamming bound of the code construction in Theorem [21 which is 

a{5*)<l-2H2{^S*). (15) 

The classical Plotkin Bound says that 

a'{6') < 1 - 26', if < 5' < i, 
a'ls') = 0, if I < 5' < 1, 

and by (fT2l) and (fT3|) . the corresponding quantum Plotkin bound of the code construction 
in Theorem [2] is 

a{6*)<l-46*, if0<6*<l- 

a{6*) =0, if i < 5* < 1. ^ ' 

The classical Elias Bound says that 



a'i5') < 1 - H,il - ^\{\ - 5')), if < 5' < i, 
a'{5') = 0, if i < 5' < 1, 
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and by (fT2l) and (fT3l) . the corresponding quantum Elias bound is 



ai6*) <l-2H^{\~ ^\{\ - 5*)), if < 5* < i; 
a{5*) = 0, if i < 5* < 1. 

The classical weaker McEliece-Rodemich-Rumsey- Welch (MRRW) bound says that 



(17) 



a\5') < H,{\ - ^5'{l - 5')) 



and by ( |T2l) and ( |T3l) . the corresponding weaker quantum MRRW bound is 



a 



(S*) < 2H,{\ - ^5*{l-5*)) - 1, if i < H,{\ - ^5*{l-5*) ) < 1; 
«(5*) = 0, if0<//2(i- V5*(l-5*))<|. 

The classical singleton bound says that 

a\6') < 1 - 5', 

and by fll2p and fll3l) . the corresponding quantum singleton bound is 

a{5*)<l-25*. (19) 
The classical Gilbert- Varshamov bound says that 

and by f[T^ and f[T^ . the corresponding quantum Gilbert- Varshamov bound is 

a{5*)>l-2H2{5*). (20) 

The above asymptotic bounds for the stabilizer code construction in Theorem [2] are depicted 
in Fig. Ei 

We next compare these asymptotic bounds with known bounds of quantum codes in the 
literature. In [H HU], the quantum Hamming bound says that for an [[n, k,d > d = 2t + 1]\ 
quantum code, 



and the asymptotic form is 



-<l--\og,3~H,{^-) 
n n n 



or 



a{6)<l-^6log,3-H2{^6). (21) 

The quantum singleton bound [HI [IE] says that for an [[n, k, d]] quantum code, 

n-k>2d-2, 
17 




Figure 2: Asymptotic coding bounds for the construction stated in Theorem [2l 



or 

a{6*) <l-26. (22) 

The Gilbert- Varshamov bound for a general quantum stabilizer codes, proved in Theorem 2 
in [11], says that an [[n, k,d = 2t + 1]] stabilizer code exists if 

->l--log2 3-if2(-), 

n n n 

or 

a{5)>l-5\og2^-H2{5). (23) 

The Gilbert- Varshamov bound for CSS codes, proved in Section V in [8], says that an 
[[n, fc, (i = 2t + 1]] CSS code exists if 

^>l-2/7,(^), 
n n 

or 

a((5) > 1-2/^2 (5)- (24) 

The above known quantum bounds in the literature are depicted in Fig. [3J It can be seen 
that the two singleton bounds flTIJl) and fl22]) for the stabilizer code construction in Theorem [2] 
and for the general quantum codes, respectively, are exactly the same. And the two Gilbert- 
Varshamov bounds (l20l) and (12^ for the stabilizer code construction in Theorem [2] and 
for CSS codes are also exactly the same. The Gilbert- Varshamov bounds (!23|) for general 
stabilizer codes is still better than the Gilbert- Varshamov bounds (l20l) for the stabilizer code 
construction in Theorem O 
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Figure 3: Known quantum bounds in the literature. 



4 Quantum Reed-Muller Codes 



In this section, we will give a family of quantum stabilizer codes from the parity-check 
matrices of Reed-Muller codes by Theorem [2l Parity-check matrices of Reed-Muller codes 
are commutative by the \u\u + v\ construction. Permutation matrices that increase the 
quantum minimum distance by Theorem O are also investigated. 



4.1 Properties of Classical Reed-Muller Codes 



Reed-Muller codes are weakly self-dual codes and have simple but good structure properties 
[16] . A Reed-Muller code with two parameters r,m is denoted by RM{r,m), < r < m. 
This code is of length 2"* and r is called its order. Consider the following (m + 1) 2''"-tuples 



1 

V2 



1111 

10 1 
11 





11111111 

10 10 10 1 
1 1 1 1 

1 1 1 1 



1111), 

10 1), 
11), 

1111). 
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Then RM{r, m) is generated by 

(degree 0) 1, 
(degree 1) Vi, ...,Vm, 
(degree 2) viV2, Vm-iVm, 

(degree r) viV2 ■ ■ ■ f^, • • • , Vm-r+iVra-r+2 ■■ - Vm, 

where the product of the fj's means the bitwise AND of the fj's and the degree means the 
number of fj's appearing in the product. There are several properties of i?M(r, m) which 
can be derived directly from its construction [16]. The dimension of RM{r,m) is 




and the minimum distance of RM{r, m) is 

d = 2™-''. 

The dual code of RM{r, m) is 

RM{m — r — l,m) 

for < r < m, where i?M(— 1, m) = {0}. Let G(^r,m) denote a generator matrix of RM(r, m). 
Reed-Muller code RM{r, m + 1) can be obtained from RM{r, m) and RM{r — 1, m) by using 
the \u\u + v\ construction. A generator matrix of RM{r, m + 1) is 

Since RM{m — r — 1, m) is the dual code of RM{r, m), a parity-check matrix of RM{r, m) 

is G(^m-r-l,m) ^ud 

G{r,m)G(^^_j.__l^^-j = O. 

For convenience, the orthogonality of Reed-Muller codes can be remarked in the following 
lemma. 

Lemma 6. For r + s < m — 1 and m > 1, G(^r,m) and G(^s,m) are orthogonal. 

Next we will consider the relation between the commutativity of G(^r,m) and the param- 
eters r, m. We first consider the case when the permutation matrix P is an identity matrix 
I in (^^. It is trivial that G(o,m) is commutative for all m > 1 since it is an all 1 vector. In 
general, we have the following lemma. 

Lemma 7. For r < [yj and m > 1, G(r,m+i) is commutative. 



Proof. By ( l25l) . we have 



G 



(r,m+l) 



G{r^m) G(^r,m) 
O G(r-l,m) 
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Since r + (r — 1) ^ m — 1, by Lemma El we have 

G{r,m)G(^j._l^^^ = O. 

Then the commutativity of G{r, m + 1) can be verified as 



O G(^r-l,m)J \Gjr,m) O 

r^T _|_ r^T r^T 

o. 



□ 



4.2 Quantum Reed-Muller Codes from Parity- Check Matrices 

Now we will present a class of stabilizer codes derived from Reed-Muller codes by Theorem 

El 

Let m > 2r and let H = G(^r,m+i) be a check matrix of a stabilizer group S. Then the 
stabilizer code C{S) is a quantum code of length n = 2"^ and dimension 

A; = 2™-^("^ + ^ 

The classical minimum distance of the parity-check matrix G(^r,m,+i) is 

^1 _ 2(m+l)-((m+l)-r-l) _ 2'-+l 

Thus the quantum error- correcting capability t of C(iS) is lower bound by 

t' = = r-' - 1, (26) 

by TheoremEl Therefore, this quantum code C{S) is able to correct at least (2^'^^ — 1) qubit 
errors provided that r > 1. The quantum minimum distance d of C{S) is lower-bounded 
by 2t' -|- 1 = 2*" — 1. On the other hand, the symplectic dual has a generator matrix 
G = G(m-r,m+i)A wliich generates the same code as the generator matrix G(^rn-r,m+i) and its 
generalized Hamming weight is 2^'. 

We have the following theorem. 

Theorem 8. The parity-check matrix of a classical Reed-Muller code RM{r, m + 1) in (l25l) 
with m > 2r and r > 1 is a check matrix of a [[2™, 2"^ — Yll=o {"^t^) ' 2'']] quantum stabilizer 
code. 
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Since the quantum error correcting capability t = 2^~^ — 1 equals to the lower bound 
t' = 2^'~^ — l by Theorem[21 this quantum code will have additional correctable error operators 
of weight > t. 

Since Reed-MuUer codes RM{r, m) with 2r + 1 ^ m are weakly self-dual codes, by 
Lemma [6l we can use them to construct CSS codes. Take Ci = RM{ri,m) with minimum 
distance 2™"^!. Then choose C2 = RM{r2, m), a subcode of Ci, with r2 < ri. The dual code of 
C2 is C2 = RM{m — r2 — l,m) with minimum distance 2"^~"^~^^~^ = 2''2+^. By CSS construc- 
tion, we obtain a quantum code with parameters [[2™, X]iir2+i (T)'— iiiin{2™~''i, 2'"2"'"^}]]. 
For the best efficiency, we take r2 + 1 = m — ri. Let r = r2 + 1. Then 2r < r2 + ri + 1 = m. 
We now construct a CSS code with parameters 

[[2",2"-2xJ ("'Y2n]. (27) 



i=0 



Comparing the dimension of a CSS code in (1271) with that of a quantum code in Theorem [HI 
both having the same length 2™ and the same minimum distance 2"^, 

p--2i:(7))-(2"-i:('"r))-(:)- 

we find that the CSS construction has a higher efficiency. However, the construction in Theo- 
rem [S] gives us additional correctable error operators of weight > t = 2^~^ — 1. In Table [H we 
list the number of additional correctable error operators for the stabilizer codes constructed 
by Theorem [2] from Reed-MuUer codes with parameters (m, r) = (5, 2), (6, 3), (7, 3). 



Table 1: A list of numbers of additional correctable error operators for the stabilizer 
codes constructed by Theorem [2] from Reed-Muller codes with parameters (m, r) = 
(5, 2), (6, 3), (7.3). 



n = 2™ 


32 


64 


128 


t 


1 


3 


3 


of additional correctable error operators 


1984 


5.99E+09 


3.87E+11 


7^ of original correctable error operators 


97 


1.14E+06 


9.29E+06 


# additional/ # original 


20.45 


5.24E+03 


4.16E+04 


m 


5 


6 


7 


r 


2 


3 


3 


dimension deficit (™) relative to the CSS code 


10 


20 


35 



On the other hand, when comparing the efficiency of the [[2"^, 2"^ - YJr}^ ("+^), 2"-i]] 
stabilizer code constructed by Theorem [2] from a Reed-Muller code with that of the [[2*", 2"^ — 
2 (T)'^*^]] code, the former code has a surplus (.™^) in dimension, while the min- 
imum distance of the former quantum code is only half of that of the latter CSS code. 
However, the former code has a lot of additional correctable error operators which will 
strengthen the error performance of the former code. 
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4.3 Permutations Which Increase the Minimum Distance 



We find that if we multiply G(i m+i) by the permutation matrix P' 



with P being 



P P. 

the permutation matrix used in [10], a stabilizer group S with a check matrix H = G(i^rn+i)P' 
will give a quantum stabilizer code C{S) with parameters [[2"^, 2™ — m — 2, 3]], which are the 
same as those constructed in [10]. For example, when m = 3, 



/O 1 0\ 

1 

1 

1 

1 

1 

1 

\i oy 

There are many other permutation matrices that will work by Theorem 10 in [15], which 
means that a column permutation on a parity-check matrix may give a stabilizer code C{S) 
with higher quantum error- correcting capability. However, it is in general hard to find such a 
permutation matrix that the commutative condition still holds after the column permutation 
of the parity-check matrix. 

We now investigate the effect of permutation matrices on generator matrices of Reed- 
MuUer codes. 



Theorem 9. Let P' 



where P is a permutation matrix such that the assumptions 



7 O' 
P P. 

dn}, (Uni), and (HI]) in Theorem hold with Ci = RM{m -r - l,m), Ca = RM{m - r,m) 
and di = 2''+^, d2 = 2^. Then for m > 2r and r > 1, the quantum stabilizer code C{S) with 
a check matrix H = G(^r,m+i)P' will have parameters [[2™, 2™ - ^^'^^ C"^^), > + 2''^^]]. In 
addition, C{S) will have additional correctable error operators if r > 3. 



Proof. By Theorem |5l the minimum distance of C{S) is at least 

min{di, [^1} = min{2^+\ [^1} = r + 2^-\ 

Note that the classical minimum distance of the parity-check matrix H = G(r,m+i)P' remains 
unchanged after a column permutation. By fl26l) . t' = 2^~^ — 1. For c? = 2^ + 2''"^, the error- 
correcting capability of the quantum code C{S) is t = [^y^J = |2^'~^ — 1. Then we have 
2t' ~ t = 2^'"^ — l>0ifr>3. Thus by Theorem [3l the quantum code C{S) will have 
additional correctable error operators of weight > t = §2''"^ — 1 if r > 3. □ 



In [19], Steane gave a class of quantum Reed-Muller codes with parameters [[2™', 2™' — 
Hl=o C^t^) 5^2^ + 2^~^]] as given by Theorem [H If there exists a permutation matrix P' 
satisfying the assumptions in Theorem [HI a stabilizer group S with a check matrix H = 
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G{r,m+i)P' , r > 1, will give a quantum stabilizer code C{S) having the same parameters as 
those in [IS] but having additional correctable error operators. 

We now give an effective permutation matrix for Gn^m)- Let 



T 



1 
1 



and 



Q 



o 



O Ir_ 



1 

1 



where n = 2™. For example, when n = 2^ = 8, 

/I 
10 



T 





10 
1 

VOOOO 0001/ 



0\ 



10 

10 





10 



/I 



1 



Q 










1 
1 









0\ 






1 

1 

1 









1 

0/ 



\0 

It can be verified that ViT = Vi+i for 1 < i < {m — 1) and VmT = vi. Thus 

P = TQ 

is a permutation matrix such that ViP = for 1 < i < {m — 1) and VmP = Vi + Vm- 

Theorem 10. The permutation matrix P = TQ is a permutation matrix such that the 
assumptions iQ, ffTOj) . and ffTTj) in Theorem [5] hold with Ci = RM{m — 2,m), C2 = RM{m — 
1, m) and di = 4, ^2 = 2 for m > 2. Then a quantum stabilizer code C{S) with a check 
matrix H = G(i,m+i)-P' will have parameters [[2*^, 2™ — m — 2, > 3]]. 



Proof. Hi = (^(i,™) and H2 = G(o,m) are parity-check matrices of Ci = RM{m — 2,m) 
and C2 = RM{m — l,m), respectively. It can be verified that HiPH-f = HiP^H^ and 
H1PH2 = O. Denote ViP = f • for convenience. It is obvious that (wj^ ■ ■ ■v^)P = f,'^ " " ""f^i, 
for any /. Since C2 = RM{m — 1, m) is generated by 

l,Vi,..., Vm, V1V2, . . . , Vm-iVm, • • • , V1V2 " " " ^m-l, . . . , V2V3 ■■■Vm, 
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C2P = RM{m — 1, m)P is generated by 

IP, ViP, VmP, {VlV2)P, ...,{ )P, . . . , {ViV2 ■ ■ ■Vm-l)P, • • • , {V2V3 ■ ■ ■Vm)P 

or 1,V2, . . .,Vm,Vi+Vm, V2V3, . . .,VmVl+Vm, {V2V3 ■ ■ ■ Vm) , ■ ■ ■, ■ " " Wm^^l) + (f 3f 4 " ' ' Vm) ■ 

It can be easily verified that C2P = C2. Similarly, we have CiP = Ci. Two generator matrices 

r " 

of Ci = RM{ni — 2,m) and C2 = RM{m — l,m) are Gi = G(m-2m) and G2 = ^ 

L Gi 

G(m-i,m)5 respectively. The row vectors of Gs are V\V2 ■ ■ ■ fm-i, • ■ •, V2V'i ■ ■ ■ Vm, the m gener- 
ators of degree m — 1. For convenience, we denote f ■ ■ ■ f ^ = Wj if j is not equal to any 
ii, 1 < / < m — 1. Then the row vectors of G3 are Wm, Wm-i, • • • , W2, Wi and the row vectors 
of G3P are WmP, Wm-iP, • • • , W2P, wiP with 

WmP = iViV2 ■ ■ ■ Vm-l)P = V2V3 ■ ■ ■ f„ = Wi, 
Wm-lP = V2V3 ■ ■ ■ Vm-lVm + V2V3 ■ ■ ■ Vm-lVl = Wi + Wm 
Wm-2P = V2 - ■■ Vm-2Vm + V2 ■ ■ ■ Vm-2VmVl = ^2 " " " ^^m-2^m + W^m-l 



W2P = V2Vi ■ ■ ■ Vm-lVmVm + V2V4, ■ ■ ■ VmVl = V2V4, ■ ■ ■ Vm + W3, 
WiP = V3V4^ ■ ■ ■ Vm-lVmVm + ^sV^ ■ ■ ■ V^Vx = VsV^ ■ ■ ■ + W2. 

If uG^P = UG3 for some binary m— tuple u = (ai, . . . , am), 

m m 
1=1 i=l 

It can be verified that the above equation holds only if ai = 02 = . . . = flm = 0, that is, 
u is a zero tuple. Therefore, uG^P 7^ nGa for any nonzero (^2 — /ci)— tuple u. By Theorem 
[9l we have a [[2"^, 2"^ — m — 2, > 3]] quantum stabilizer code C{S) with a check matrix 
H = G(i,^+i)P' for m ^ 2. □ 



There are other permutation matrices which will work by similar proofs. 

For a general Reed-MuUer code RM{r,m) with minimum distance d' = 2^""^, m > 2r, 
we have the following conjecture. 

Conjecture 11. Either P = T ot P = TQ is a permutation matrix satisfying the assump- 
tions in Theorem [9] for all r, m with m > 2r. □ 



5 Quantum Cyclic Codes 
5.1 Quantum Circulant Codes 

Motivated by the five-qubit code in Eq. ([2]), we use Construction II in Subsection 13.21 with a 
check matrix H = [Hx\Hz] = [G1IG2], where Gi,G2 are generator matrices of two classical 
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cyclic codes €1,62, respectively. This method is called the circulant construction and the 
generated quantum codes are called quantum circulant codes. 

We arbitrarily choose two binary polynomials gi{x),g2{x) as generator polynomials of 
two classical cyclic codes €1,62, respectively. Then they cyclicly generate the matrices 
Hx = Gi,Hz = G2, respectively. If H is commutative, we then justify the rank of H 
by transforming H into the row-reduced-echelon form. The minimum distances of these 
quantum codes are determined by computer search and for small n, many good quantum 
circulant codes are found by this method and the quantum circulant codes with parameters 
achieving the upper bound in Table III in [15j are listed in Table [21 It is conjectured that 

Table 2: Some extremal quantum circulant codes. 



n\d\k 





1 


2 
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5 
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3 
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4 
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2 
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3 






2 


2 


2 
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3 


3 


3 






2 






10 


4 


4 








2 


2 
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12 


6 
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4 
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13 
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5 
















14 
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4 


4 
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15 
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5 


5 


5 


4 


4 
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3 
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16 


6 


6 








4 








17 


7 


7 
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18 






6 




5 










19 




7 

















most of the best quantum codes can be obtained by the circulant construction. However, the 
minimum distance is hard to decide for large n for the extremely high computing complexity. 

5.2 Quantum Quadrature- Residue Codes 

Another interesting fact of the five-qubit code is that the binary bit patterns of the generator 
polynomials gi{x),g2{x) are the indicator vectors of the quadratic-residues mod 5 and the 
non- residues, respectively. Inspired by this fact, a suggestion for any prime p of the form 
8m + 5 is given in [TT] without a proof though there is no classical binary quadratic-residue 
code for theses p's. In Subsubsection 15.2. 1^ we first introduce the basics of quadratic residues 
(please refer to [TB] for a detailed discussion) and then give several quantum quadratic-residue 
codes by the CSS construction for prime numbers of the form p = 8m ± 1. We then use 
the circulant construction with the two indicator vectors of the quadratic-residues and the 
non-residues mod p for prime numbers p of the form p = Aj + 1 in Subsubsection I5.2.2[ The 
minimum distance of the generated quantum codes will be calculated by computer search. 
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5.2.1 Binary Quadrature-Residue Codes and Their Applications 



Let p be an odd prime. For a nonzero j,l<j<p — 1, ifj — P mod p for some I, we 
say that j is a quadratic-residue mod p, otherwise j is called a non-residue. In particular, 
is neither a residue nor a non- residue. There arc — 1) quadratic-residues and |(p — 1) 
non-residues mod p. Let Q denote the set of quadratic residues mod p and denote the set 
of non-residues. If p is a primitive element of the field GF{p), then />* e Q if and only if j is 
even. Thus Q is a cyclic group generated by p^. It can be verified that Q is a disjoint union 
of cyclotomic cosets mod p. Let 

q{x) = Yl{x - a''), n{x) = Yl{x -a""), 

reQ neN 

where « is a primitive p-th root of unity. Then x^—l = {x—l)q{x)n{x). The quadratic-residue 
(QR) codes Q, Q,J\f,J\fa.Te cyclic codes with generator polynomials q{x), {x—l)q{x),n{x), {x— 
l)n{x), respectively. The dimensions of Q,J\f are both ^{p+1) and the dimensions of Q,Af 
are both |(p - 1). It is clear that Q,Af ^J^. 

A binary polynomial E{x) in the quotient ring Rp = GF{2) [x]/ {x^ — 1) is an idempotent 
if E{x) = E{xY = E{x^) mod x^ — 1. A binary cyclic code C — {g{x)) of length p contains a 
unique idempotent E{x) such that C — {E{x)). Of course, E{x) — p{x)g{x) for certain p{x) 

Consider a prime p such that 2 e Q. li p = Aj — 1, the a can be chosen such that the 
idempotents of Q,Q.,N,M are 

E,{x) = J2 ^^ F,{x) = 1 + ^ x", E^{x) = J2 ^n(^) = 1 + E 

reQ neN neN reQ 

In this case, = Q and Af^ = Af. If p = 4j -|- 1, the a can be chosen such that the 
idempotents of Q, Q,Af,Af are 

E,(x) = 1 + E F,{x) = ^ x", En(x)^l + J2 Fn(x) = a;^ 

reQ neN neN reQ 

In this case, ^ Af and Af^ = Q. 

If Fq[x) = Yl^=o fi^^j then a generator matrix of Q is 



fo 


/l ■ 


• fp-l 


fp-l 


fo . 


■ fp-2 


/l 


/2 . 


■ fo 



and a generator for Q is 

G 

Similar results can be obtained for Af and Af. 
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It is obvious that the above dual pairs can be used in the CSS construction. Taking 
Ci = Q, C2 = Q or Ci = A/", C2 = A/" if p = 4j - 1, and Ci = Q, C2 = TV" or Ci = A/", C2 = Q if 
p = 4j + 1, we can obtain a [[p, 1, d]] quantum code for certain d in both cases. By Theorem 
1 and Theorem 23 in Ch.l6 in [16], we have the following theorem similar to Theorem 40 
and Theorem 41 in [20] but of binary case. 

Theorem 12. li p = 8m ± 1, then there exists a [[p, l,(i]] CSS code with d > y^. If 
p = 4j — 1, we can strengthen the bound to d"^ — d + 1 > p. 

By collecting the minimum distance of various classical binary quadrature-residue codes 
in [211 ESI [23] , we obtain Table El 

Table 3: Parameters of several [[p, 1, d]] CSS codes from classical binary quadrature-residue 
codes. 



p 


7 


17 


23 


31 


41 


47 


71 


d 


3 


5 


7 


7 


9 


11 


11 


p 


73 


79 


89 


97 


103 


113 


137 


d 


13 


15 


17 


15 


19 


15 


21 



Quadrature-residue codes can be extended by adding an overall parity-check bit so that 
the extended quadrature-residue codes Q and A/ have the following relations: 

Q-^ = Q and A^^ = A^, if p = 4j - 1 

and 

Q^=Af, ifp = 4j + l. 

Similarly, the extended quadrature-residue codes can be used in the CSS construction and 
we have the following theorem by Theorem 8 in Ch.l6 in [T6]. 

Theorem 13. If p = 8m ± 1, then there exists a [[p -|- 1, 0, d]] CSS code for certain d. If 
p = 4j — 1, d = or 3 mod 4. If p = 4j -|- 1, d is even. 

With Table 1(a) in [24j, we have Table H 

The parameters in Table [3] are related to those in Table H] by Theorem 6 in [15] in spite 
of the fact that the entries are fewer in Table [31 

5.2.2 Quadrature Residues Related Quantum Circulant Codes 

Let gi{x) = YnlQaiX' = XlreQ^'" 92{x) = J^jZlbjX^ = Y.n&N^"' certain prime 
number p. These two binary polynomials are corresponding to the indicator vectors of the 
quadratic-residues and the non-residues, respectively. For example, when p = 13, 

g^[x) = X + + x"^ + + + x^^ g2{x) = x^ + x^ + x^ + x'^ + x^ + x^^ 
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Table 4: Parameters of several [\p + 1, 0, d\] CSS codes from classical binary extended 
quadrature-residue codes. 



p 


8 


18 


24 


32 


42 


48 


72 


d 


4 


6 


8 


8 


10 


12 


12 


p 


74 


80 


90 


98 


104 


114 


128 


d 


14 


16 


18 


16 


20 


16 


20 


p 


138 


152 


168 


192 


194 


200 




d 


22 


20 


24 


28 


28 


32 





and the two corresponding indicator vectors are 

0101100001101 , 0010011110010, 
respectively. We begin with the following lemma to discuss our method for p — Aj + 1. 
Lemma 14. For p — Aj + 1, the matrices 



and 



are symmetric. 



Hp 



ap-2 



ai 


a2 . 


ao 


bo 


bi . 


■ bp_i 


'p-i 


bo . 


■ bp_2 


bi 


b2 . 


. bo 



Proof. We have Uj = 1 if j e Q and aj = 0, else, by the choice of gi{x). For p = Aj + 1, 
— 1 e Q and we have Uj — a^j — Up^j. The element in the i-th row and the j-th column of 



is H 



X)ij 



mod p- Since (yHx)ij (^j—i mod p ^—(j—i) mod p 



^i—j mod p {,Hx)jii 
□ 



the matrix Hx is symmetric. Similar reason holds for Hz- 

Theorem 15. For p = 4j ' + 1, the matrix H = [Hx\Hz] is commutative and has rank p—1. 
Then there is a [[p, 1, d]] quantum stabilizer code for a certain d. 



Proof. Let Hy = Hx + Hz. Equivalently, we consider a check matrix [Hy\Hz] instead of 
[Hx\Hz]. Note that the rank of the check matrix remains to be determined, 
the p X p all-1 matrix. It is obvious that 



Let Jp denote 



H. 



1 1 

1 1 
1 1 

111 



Jp I p. 
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Denote the i-th row vectors of Hy and Hz by ai and respectively. Since and /3j are 
the indicator vectors both with i right cychc shifts, we have ■ j3i = 0. Note that a, + 
aj = . . . 010 . . . 010 ... with two I's at the i-th and the j-th positions. The commutative 
condition can be checked for any two rows as follows: 

ai ■ (3j + aj ■ f3i = (a, + aj) ■ {(3i + (3j) 

= (0 ... 010 ... 010 ... 0) ■ + Pj) 

= {Hz)i,i + {Hz)j,i + {Hz)i,j + {Hz)j,j . 

Since Hz is symmetric by Lemma [H] and {Hz)i,i = for all i by the construction, the 
commutative condition holds. For Hy = Jp — Ip, the last row vector of Hy can be obtained 
from the summation of all the other row vectors. Thus the rank of Hy is p — 1. Similarly, 
the ranks of Hx and Hz are both p—1. Hence the rank of the check matrix H is p — 1. Thus 
the quantum stabilizer code with a check matrix H has parameters [[p, l,(i]] for a certain 
d. □ 



In [TT], the quantum codes of p = 8m + 5 is a special case of above theorem. It remained 
to determine the minimum distances of these codes. Parameters of several quantum codes 
from Theorem [15] are given in Table [51 The minimum distance of these codes is determined 
by computer search. Note that the minimum distance of the quantum codes in Table [S] with 



Table 5: Parameters of several [[p, l,d]] quantum codes from Theorem [T5l 



p 


5 


13 


17 


29 


d 


3 


5 


5 


11 



p = 5, 13, 29 achieves the upper bound in [15]. However, the minimum distance of a generic 
quantum code from Theorem [15] is not found for the extremely high computing complexity. 



5.2.3 A Construction for Quantum Codes w'lih k = 1 



Inspired from the proof of the quantum quadratic residue codes, we give a construction of 



[[n, 1]] quantum stabilizer codes in this subsubsection. Let g = (oq, cti, . 
of length n, odd or even, with = and a, = a„_j = a_j mod n for i 
n X n matrices Hx, Hz for convenience of explanation. Let 



, ttn-i) be a vector 
1 to We use 



H 



X 



1 
1 








1 
1 



1 1 





In-l 



1 





with {H^ 



X, 



{Hx)i n-i = I ioT < i < n — 2 and {Hx),i j = otherwise. 
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Let Hz 



+1) mod n ~l~ ^(i—j) mod 


„ for 0<i,j 


< n — 


l.That 


is 


Oi + ao 


O2 + C^n— 1 


+ an-2 


• • • a„ 


-1 + 02 


Oi 


Oi + «! 


0.2 + «0 




• • • '^n 


-1 + 03 


02 


ai + 02 


02 + 0-1 


O3 + Oo 


• • • (^n 


-1 + 04 


03 


(2i ~r (Xr) 9 


do ~r Gf"n 


O3 ^ 0„_4 


■ ■ ■ o„ 


-1 Oq 


On-l 


/I ^ 1 rt -, 

U-l U-n-l 


(^2 ~\~ 0,n—2 


03 + 0"n-3 


• • • a„ 


-1 + Oi 


Oo 


ai 


(^2 + 0,n-l 


n n A- n n 


• ■ ■ 


-1 "'2 


LLl 


ai + ai 


02 


(^3 + C^n-l 


• • • 


-1 + ^3 


O2 


ai + 02 


O2 + Oi 


O3 


■ ■ ■ 


-1 + 04 


O3 


oi + a„_2 


O2 + 0„_3 


^3 + 0,n-A 






On-l 




















Then a check matrix H = [Hx\Hz] has rank {n — 1) and the commutative condition for H 
can be justified as follows: 



(0...010...010...0) ■ iPi + Pj] 



H 



Z{i,i) 



+ Hz{ij) + H 



O'i+l + (o-i-i mod n + + ('^j-i 

0. 



mod n + fli+l) + 



We found that the quantum quadratic-residue code for p = 8j + 5 in the last subsubsection 
can be constructed in this way with g being the indicator vector of the quadratic residues. 
Some quantum codes achieving the upper bound in [15] can be constructed similarly. For 
example, when n = 17, each of the following two vectors 

0110100110010110, 0100011111100010 

(or their complementary vectors and no others) together with oq = gives a [[17,1,7]] 
code that achieves the upper bound in [TJ]. A [[17,1,7]] quantum stabilizer code can be 
constructed by quantum BCH codes in [251 ES]- However, the above two vectors are found 
by computer search. It is difficult to determine a vector g and the minimum distance of the 
resulted quantum code efficiently. 



6 Conclusion 



In this paper, a simple stabilizer code construction was proposed based on syndrome assign- 
ment by classical parity-check matrices. The construction of quantum stabilizer codes can 
then be converted to the construction of classical binary linear block codes with commutative 
parity-check matrices. The asymptotic coding performance of this construction was shown 
to be promisingly comparable to that of the CSS construction. 
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Permutation matrices may help transform non-commutative parity-check matrices to 
commutative ones and/or increase the minimum distance of the constructed quantum codes. 
However, for a given parity-check matrix H , it remains open to find an effective permutation 
matrix P such that HP is commutative and/or corresponds to a code with greater minimum 
distance. 

We have constructed a family of stabilizer codes from classical binary Reed-Muller codes 
with performance comparable to that of the CSS construction. We have also investigated 
sufficient conditions for permutation matrices to be able to increase the minimum distance of 
our constructed quantum Reed-Muller codes by half. We have also proposed a specific kind of 
effective permutations and showed that they meet the sufficient conditions for stabilizer codes 
constructed from the RM{l,m) Reed-Muller codes. A conjecture of effective permutation 
matrices for general r, m with m > 2r remains to be proved. It is believed that permutation 
of the columns of a parity-check matrix will play an important role in the construction of 
quantum stabilizer codes from classical parity-check matrices. 

The quantum quadratic-residue codes are codes with large quantum minimum distance. 
However their quantum minimum distance is hard to determine as in the classical case. 
Perhaps techniques in [23] can be applied to the quantum case. How to determine the 
minimum distance of a long quantum code will be a key to find good codes. 



References 

[1] p. W. Shor, "Rapid communications," Physical Review A, vol. 52, no. 4, pp. 2493-2496, 
1995. 

[2] A. M. Steane, "Error correcting codes in quantum theory," Physical Review Letter, 
vol. 77, no. 5, pp. 793-797, 1996. 

[3] , "Simple quantum error-correcting codes," Physical Review A, vol. 54, no. 6, pp. 

4741-4751, 1996. 

[4] A. Ekert and C. Macchiavello, "Quantum error-correction for communication," Physical 
Review Letter, vol. 77, no. 12, pp. 2585-2588, 1996. 

[5] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. Wootters, "Mixed state 
entanglement and quantum error correction," Physical Review A, vol. 54, no. 5, pp. 
3824-3851, 1996. [Online]. Available: | http://arxiv.org/abs/quant- ph/9604024l 

[6] E. Knill and R. Lafiamme, "A theory of quantum error- correcting codes," Physical 
Review A, vol. 55, no. 2, pp. 900-911, 1997. 

[7] R. Lafiamme, C. Miquel, J. P. Paz, and W. H. Zurek, "Perfect quantum error correcting 
codes," Physical Review Letter, vol. 77, no. 1, pp. 198-201, 1996. 



32 



[8] A. R. Calderbank and P. W. Shor, "Good quantum error- correcting codes exist," 
Physical Review A, vol. 54, no. 2, pp. 1098-1105, 1996. [Online]. Available: 
|littp : /7arxi V. org/abs/ quant- ph/ 95 1 2032] 

[9] A. M. Steane, "Multiple particle interference and quantum error correction," 
Proa R. Soc. London A, vol. 452, pp. 2551-2576, 1996. [Online]. Available: 
[http: / 7arxiv.org/abs/quant-pli/9601029| 

[10] D. Gottesman, "Class of quantum error-correcting codes saturating the quantum 
hamming bound," Physical Review A, vol. 54, no. 3, pp. 1862-1868, 1996. [Online]. 
Available: [http: //arxiv.org/abs/quant-ph /9604038| 

[11] A. R. Calderbank, E. M. Rains, P. W. Shor, and N. J. A. Sloane, "Quantum error 
correction and orthogonal geometry," Physical Review Letter, vol. 78, no. 3, pp. 
405-408, 1997. [Online]. Available: ' http://arxiv.org/abs/quant-p h/9605005 

[12] D. Gottesman, "Stabilizer codes and quantum error correction," Ph.D. dissertation, 
California Institute of Technology, Pasadena, CA, 1997. [Online]. Available: 
[http : / / arxi V. org / abs / quant- ph/9705052 

[13] A. M. Steane, "Enlargement of Calderbank-Shor-Steane quantum codes," IEEE 
Trans. Inform. Theory, vol. 45, no. 7, pp. 2492-2495, 1999. [Online]. Available: 
http: / / arxiv.org/abs/quant-ph/980206T] 

[14] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information. 
Cambridge, UK: Cambirdge University Pressr, 2000. 

[15] A. R. Calderbank, E. M. Rains, P. W. Shor, and N. J. A. Sloane, "Quantum error 
correction via codes over GF{A)," IEEE Trans. Inform. Theory, vol. 44, no. 4, pp. 
1369-1387, 1998. [Online]. Available: ]http://arxiv.org/abs/quant-ph /9608006] 

[16] F. J. MacWilliams and N. J. A. Sloane, The Theory of Error- Correcting Codes. Ams- 
terdam, The Netherlands: North-Holland, 1977. 

[17] J. H. van Lint, Introduction to Coding Theory. New York: Springer- Verlag, 1982. 

[18] J. Preskill, "Advanced mathematical methods of physics - quantum computation and 
information," Physics, vol. 229, 1998. 

[19] A. M. Steane, "Quantum Reed-Ruller codes," IEEE Trans. Inform. 
Theory, vol. 45, no. 5, pp. 1701-1703, 1999. [Online]. Available: 
http:/ /arxiv.org/abs/arXiv:quant-ph/9608026| 

[20] A. Ketkar, A. Klappenecker, S. Kumar, and P. K. Sarvepalli, "Nonbinary stabilizer 
codes over finite fields," IEEE Trans. Inform. Theory, vol. 52, no. 11, pp. 4892-4914, 
2006. 

[21] D. Coppersmith and G. Seroussi, "On the minimum distance of some quadratic residue 
codes," IEEE Trans. Inform. Theory, vol. 30, no. 2, pp. 407-411, 1984. 



33 



[22] N. Boston, "The minimum distance of the binary [137, 69] quadratic residue code," 
IEEE Trans. Inform. Theory, vol. 45, no. 1, p. 282, 1999. 

[23] T.-K. Truong, Y. Chang, and C.-D. Lee, "The weight distributions of some binary 
quadratic residue codes," IEEE Trans. Inform. Theory, vol. 51, no. 5, pp. 1776-1782, 
2005. 

[24] M. Grassl, "On the minimum distance of some quadratic- residue codes," in Proc. IEEE 
ISIT '2000, p. 253. 

[25] M. Grassl and T. Beth, "Quantum BCH codes," in Proc. X. International Symposium 
on Theoretical Electrical Engineering, Magdeburg, 1999, pp. 207-212. [Online]. 
Available: http: / / arxiv.org/ abs/quant-ph/99 10060 

[26] S. A. Aly, A. Klappenecker, and P. K. Sarvepalli, "Primitive quantum BCH codes 
over finit e fields," in Proc. IEEE ISIT '2006, pp. 1105-1108. [Online]. Available: 
http: / / arxiv.org/ abs/quant-ph/0501126, 



34 



